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Abstract 

The following problems are discussed. 1. Pulsars and close bina- 
ries. 2. Hulse- Taylor pulsar. 3. Disrupted pulsar pairs. 4. RP statis- 
tics. 5. Enhanced evaporation: formation of single RP. 6. General 
relativity effects: NS-I-NS. 7. A Double pulsar system. 8. Checking 
general relativity. 9. Variability of the gravitational constant. 10. 
Space Watch. 



PACS numbers: 97.60.Gb; 95.30.Sf. 
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1 Introduction 

Discovery of pulsars was announced by Hewish, Bell et al.[l]. Until 1973 all 
known pulsars (more than 100) had been single, while more than half massive 
stars (predecessors of pulsars) are in binaries. Possible explanations of this 
fact had been suggested: all pairs are disrupted during explosion, or there is 
no possibility to form a radiopulsar in pairs [2]. 



*Space Research Institute RAN. Profsoyuznaya 84/32, Moscow 117997, Russia; and 
JINR, BLTP, Dubna, Russia 



1 



2 Pulsars and close binaries 



X-ray satellite UHURU was launched in 1970, and soon after that X-ray 

pulsars in binaries had been discovered. One of the best studied is X ray 
pulsar Her X-1. It has pulsation period Pp ~ 1.24 sec, orbital period Port ~ 
1.7 days, neutron star mass about 1.4 Mq, and optical star with a mass about 
2Mq, see e.g [3]. 

It was shown in [4], that this system should give birth to the binary 
radiopulsar, by following reasons. 1. After 100 million years the optical star 
will become a white dwarf, mass transfer will be finished, and the system will 
be transparent to radio emission. 2. X ray pulsar is accelerating its rotation 
due to accretion, so after the birth of the white dwarf, the neurton star will 
rotate rapidly, with Pp about 100 ms. It was suggested in [4], that binary 
radiopulsars had not been not found until 1973, because the magnetic field 
of the neutron star is decreasing about 100 times during the accretion, so 
binary radiopulsars are very faint objects. Pulsar luminosity L ~ B^/P'^, so 
at small B, luminosity L is low even at the rapid rotation. Magnetic field is 
decreasing due to screening by the infalling plasma. 



3 Hulse- Taylor pulsar 

Hulse and Taylor had discovered [5] the first binary radiopulsar with a period 
Pp = 0.05903 s, orbital period Porb ~ 7.75 h, orbit eccentricity e = 0.615 (see 
Fig.l). The properties of the first binary pulsar coincide with our predictions: 
rapid rotation together with a small magnetic field [6], [7] 



(1) 




The average magnetic fields of single radiopulsars is about 10^^ Gauss. 
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Figure 1: Velocity curve for the first binary pulsar. Points represent mea- 
surements of the pulsar period distributed over parts of 10 different orbital 
periods. The curve corresponds to the pulsar parameters from the text. 

4 Disrupted pulsar pairs 

Suggestions for separation of pulsars at birth was done in [4] where 10 pairs 
of single pulsars with a possible origin from one pair had been listed (see 
Table 1). The same idea was recently considered in [8], who measured pulsar 
proper motions. They suggested a common origin of the pair B2020-I-28 and 
B2021+51 (see Fig.2). It is interesting, that one of this pulsars (2020+28) had 
been suggested in [4] for the common origin with another pulsar (2016+28, 
line 8 in Table 1). The last pulsar is much closer to B2020+28 on the sky, 
and both are situated at the same distance [9]. The pulsar B2021+51 is much 
farther on the plane, and almost two times by the distance, but velocities 
intersect in Cygnus OB association. 

5 Recycled pulsar (RP) statistics 

Recycled pulsars is a separate class of neutron stars, containing more than 180 
objects. All passed the stage of accreting pulsars, accelerating the rotation 
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Table 1: Single pulsars of possible common origin, from [4]. 



No 


PSR 


R. A. 


Decl 


1 T T 


1 T T 




DM ^ 


^ day 


p 


1 


P0943+10 


9''43™20" 


10°05'33" 


225.4 


43.2 


1.098 


15 


- 


- 




P0950+08 


9 50 31 


08 09 43 


228.9 


43.7 


0.253 


3 


0.0198 


3.5-10^ 


2 


P0809+74 


8 09 03 


74 38 10 


140 


31.6 


1.30 


6 


0.014 


2.5-10^ 




P0904+77 


9 04 


77 40 


135.3 


33.7 


1.58 








3 


PI 700- 18 


17 00 56 


-18 


4.0 


14.0 


0.802 


<40 


0.154 


6.9-106 




PI 706- 16 


17 06 33 


-16 37 21 


5.8 


13.7 


0.653 


25 


0.55 


3.3-106 


4 


P0329+54 


3 29 11 


54 24 38 


145.0 


-1.2 


0.714 


27 


0.177 


1.1-10^ 




P0355+54 


3 55 00 


54 13 


148.1 


0.9 


156 


55 






5 


P0525+21 


5 25 45 


21 55 32 


183.8 


-6.9 


3.7455 


51 


3.452 


3-106 




P0531+21 


5 31 31 


21 56 55 


184.6 


-5.8 


0.03313 


57 


36.526 


2.5-10^ 


6 


P1426-66 


14 26 34 


-66 09 94 


312.3 


-6.3 


0.787 


60 








P1449-65 


14 49 22 


-65 


315.3 


-5.3 


0.180 


90 






7 


P1845-01 


18 45 


-01 27 


31.3 


0.2 


0.660 


90 








P1845-04 


18 45 10 


-04 05 32 


28.9 


-1.0 


0.598 


142 






8 


P2016+28 


20 16 00 


28 30 31 


68.1 


-4.0 


0.558 


14.16 


0.01 


1.2-10^ 




P2020+28 


20 20 33 


28 44 30 


68.9 


-4.7 


0.343 








9 


P2111+46 


21 11 41 


46 36 


89.1 


-1.2 


1.015 


141.4 








P2154+40 


21 54 56 


40 00 


90.5 


-11.5 


1.525 


110 






10 


P0611+22 


06 11 10 


22 35 


188.7 


2.4 


0.335 


99 








P0540+23 


05 40 10 


23 30 


184.4 


-3.3 


0.246 


72 
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and decreasing the magnetic field. So we have ordinary pulsars with Pp = 
0.033 - 8 s, B = 10^^ — 10^^ Gs; and recycled pulsars with Pp = 1.5 - 50 ms, 
B — 10^ — 10^° Gs. These two groups are distinctly visible on the P — P 
diagram (Fig.3). 

6 Enhanced evaporation: formation of single 
RP 

Most recycled pulsars consist of pairs NS + WD, and single NS, about 180 
objects in total. Only 6 NS + NS pairs exist, like Hulse- Taylor pulsar. Single 
recycled pulsars (small Pp and B) have passed the stage of X-ray pulsar in 
binary, and later had lost the WD companion. Besides, NS + NS are in the 
Galaxy disk, and NS + WD, and single are in the galactic bulge, and in the 
globular clusters. Globular clusters contain 0.001 of the mass of the Galaxy, 
and about 1/2 of the recycled pulsars. Formation of the binary in GC may 
happen by tidal capture or triple collision. Disruption of the binary RP also 
happens by collisions with GC stars: most single RP are situated in GC. 
In the most dense GC we have [10]: 31 RP in Terzian 5 (11 BP, 10 +10? 
single); 22 RP in 47 Tuc (14 BP, 7+1? single). In total: 80 RP in Gal 
(15 single), 108 RP in GC (40+15? single). 

Simple disruption of the pair by collisions with field stars does not work, 
because hard pairs become even harder by collisions. When pair is in a 
state of a disk accretion, with WD fiUing its Roche lobe (Fig. 4) situation is 
opposite, and also hard pairs become softer during collisions due to enhanced 
mass transfer, and finally pair is disrupted. This process named "enhanced 
evaporation" , first considered in [1 1] , could explain formation of single RP in 
GC, and their appearance in bulge may be connected with a full evaporation 
of GC (see also [12]). 

7 General relativity effects: NS+NS 

In the NS+NS binary tidal effects are negligibly small, so several post- 
keplerian (PK) GR effects had been measured from pulsar timing. They 
include [13] parameter ui, determining the rate of relativistic apse line mo- 
tion, parameter 7, representing the amplitude of the signal time delay due to 
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Figure 2: The 3-dimensional pulsar motion through the Galactic potential for 
B2020+28 and B2021+51. The dashed circle represents the Cygnus super- 
bubble, while the labelled solid eUipses are the Cygnus OB associations with 
positions and extents. The extent of OB 2 is unknown and only the centre of 
the association is indicated. The thick solid lines indicate the pulsar paths, 
with the origin denoted by the starred symbol and the arrows pointing in the 
direction of motion. The current positions are indicated by the solid dots. 
The elliptical contours around the pulsars' origin in these panels indicate the 
1, 2 and 3a levels of the likehhood solution for the birth location, from [8] 
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variable gravitational redshift and time dilation (quadratic Doppler effect), 
when the pulsar moves in an elliptical orbit. Emission of gravitational waves 
results in the loss of the orbital angular momentum and decreases the orbital 
period Porb- The parameters r and s determine the time delay due to the 
Shapiro effect and are related to the companion's gravitational field. 

(p x —5/3 -1 
^) Y^A^A + Msf/^ (2) 

_^2/3 fPorbV^'' Mb{Ma + 2Mb) , . 



P 1927r 5/3 (Porb\ (l + je' + je') MaMb 

5 ® Utt; (l-e2)V2 (M^ + Mb)V3' 

t^TqMb, = GMq/c^ = 4.925490947 //s, (5) 

NS + NS RPs are the best laboratories for checking of General Relativity. 
1913+16 timing had shown (indirectly) the existence of gravitational waves, 
Nobel Prize of Hulse and Taylor (1993) 



8 A double pulsar system 

A binary system, containing two pulsars with periods 23-ms (J0737-3039A), 
and 2.8-sec (B) was discovered in 2004 year [14]. This highly-relativistic 
double-neutron-star system allows unprecedented tests of fundamental grav- 
itational physics. It was observed a short eclipse of J0737-3039A by J0737- 
3039B, and orbital modulation of the flux density and pulse shape of J0737- 
3039B, probably due to the influence of J0737-3039As energy flux upon its 
magnetosphere. These effects will allow to probe magneto-ionic properties 
of a pulsar magnetosphere. The observational properties of this system are 
hsted in Table.2, and artistic view is give in Fig. 5. 



7 



-12 



-14 



-16 



O 



3-lE 



-20 



-22 



10^ y . 



10^ y^. 



lO'' y^ 









-a 



■-J0l° G 



10° G 



10 100 1000 

Period (milliseconds) 



Figure 3: The location of pulsars on the P — P diagram (period-period 
derivative). Pulsars in binary systems with low-eccentricity orbits are en- 
circled, and in high-eccentricity orbits are marked with ellipses. Stars show 
pulsars suspected to be connected with supernova remnants [10]. 
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Figure 4: Artistic vision of a disc accretion to the neutron star in X-ray 
sources (http:/ /antwrp. gsfc.nasa.gov/apod/ap991219. html) 

9 Checking general relativity 

Five equations in the system (2)-(6) contain two unknown masses. For cor- 
rect gravitational theory all curves on the plane {Ma, Mb) intersect in the 
same point. Three curves available for the pulsar 1913+16 (Fig. 6), and six 
curves for the double pulsar system (Fig. 7, the additional line R follows from 
observations of both pulsars) intersect in one point inside the error limits, 
indicating the correctness of GR on the level ~0.4% for 1913-1-16, after about 
20 years of observations, and on the level ~0.1% after one year of observations 
of the double system. 

10 Variability of the gravitational constant. 

Variations of the gravitational constant G have been measured by different 
methods [15]. The following restrictions are obtained from investigations of 
stellar and planetary orbits. Using of the Viking Lander ranging data gave 



9 



limits G/ G = (2±4) • 10 ^^yr ^. The combination of Mariner 10 and Mercury 
and Venus ranging gave [16] 

G/G = (0 ± 2) • 10-^V"\ (7) 

Lunar Laser Ranging experiments (LLR) are lasting for many years, and 
improve the estimations of G/G from G/G < 6 • lO^^^yr^^ in [17], up to 
G/G [yr~^]= (6 ± 8) • 10^^'^ in [18]. Incompletely modeled solid Earth tides, 
ocean loading or geocenter motion, and uncertainties in values of fixed model 
parameters have to be considered in those estimations. 

Recycled pulsars give now the best available timing precision. The ob- 
served changes in the binary period of PSR B1913+16 coincide within error 
bars with GR prediction of emission of gravitational waves. 




Figure 5: Artistic picture of the double pulsar system, from [19]. 

Measurements of the variations of the binary period give restrictions to the 
variation of the gravitational constant. Change of G has influence on the 
orbital motion, and only emission of the gravitational waves compete here 
with the variation of G. If we accept the correctness of GR, than we can 
use residuals (error box) of the measurements of decay of the binary period 
Pb for the estimation of variations of G [21]. It was obtained in [22], using 
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12 3 
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Figure 6: Observational constraints on the component masses in the binary 
pulsar PSR 1913 + 16. The solid curves correspond to Eqns (2)-(4) with 
measured values of 7, and Porb- The intersection of these curves at one 
point (within an experimental uncertainty of about 0.35% in Porb) proves 
the existence of gravitational waves. The dashed lines correspond to the 
predicted values of the parameters r and s. These values can be measured 
by future accumulation of observational data [13]. 

timing data from [23] 

For the pulsar PSR B1913+16 the error budget for the orbital period deriva- 
tive, in comparison with GR prediction, is given in Table 3 [23]. Supposing 
that all deviations from GR are connected with G variation, we obtain the 
upper limit for these variation as follows [24]: ^ = (4.3 ± 4.9) ■ 10^^^ yr^^. 
These estimations slightly improve the values [22]. This improvement nar- 
rows on 40% the boundary in the negative region, most important from the 
theoretical point of view. Combination the last result with [16] permits to 
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Table 2: Observed and derived parameters of PSRs J0737— 3039A and B. 
Standard errors are given in parentheses after the values and are in units 
of the least significant digit (s). The distance is estimated from the disper- 
sion measure and a model for the interstellar free electron distribution, from 
[14],[20]. 



Pulsar 

Pulse period P (ms) 
Period derivative P 
Right ascension a (J2000) 
Declination 6 (J2000) 
Orbital period Pb (day) 
Eccentricity e 

Advance of periastron di (deg yr~^) 

Projected semi-major axis x = asini/c (sec) 

Gravitational redshift parameter 7 (ms) 

Shapiro delay parameter s = sin i 

Shapiro delay parameter r (/xs) 

Orbital decay A (lO^i^) 

Mass ratio R = Ma/ Mb 

Characteristic age t (My) 

Surface magnetic field strength B (Gauss) 

Spin-down luminosity E (erg/s) 

Distance (kpc) 



PSR J0737-3039A PSR J0737-3039B 
22.699378556138(2) 2773.4607474(4) 
1.7596(2) X 10-18 0.88(13) x IQ-^^ 
07'>37'"51". 24795(2) 
-30°39'40".7247(6) 
0.1022515628(2) 
0.087778(2) 
16.900(2) 
1.415032(2) 1.513(4) 

0.39(2) 
0.9995(4) 

6.2(6) 
-1.20(8) 

1.071(1) 

210 50 
6.3 X 10^ 1.6 X 10^2 

5800 X 10^° 1.6 X 10^° 

~0.6 
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Figure 7: The observational constraints upon the masses mA and mB in the 
double pulsar system J0737-3039.. The colored regions are those which are 
excluded by the Keplerian mass functions of the two pulsars. Further con- 
straints are shown as pairs of lines enclosing permitted regions as predicted 
by general relativity: (a) the measurement of the advance of periastron, giv- 
ing the total mass mA+mB = 2:588 pm 0:001M (Sun) (dashed line); (b) the 
measurement of R = mA/mB = 1:071 ± 0:001 (solid line); (c) the measure- 
ment of the gravitational redshift/time dilation parameter (dot-dash line); 
(d) the measurement of Shapiro parameter r (dot-dot-dot-dash line) and (e) 
Shapiro parameter s (dotted line). Inset is an enlarged view of the small 
square encompassing the intersection of these constraints, from [20]. 

narrow the range of G variations, which is now may be situated in the limits 
[24] 
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Table 3: Error budget for the orbital period derivative, in comparison with 
the general relativistic prediction. 





Parameter 


(10-12) 


Observed value 


pobs 


-2.4225±0.0056 


Galactic contribution 


pgal 


-0.0124±0.0064 


Intrinsic orbital period decay 


pobs pgo-l 
^b ~ ^b 


-2.4101±0.0085 


General relativistic prediction 


PGR 


-2.4025±0.0001 



G/G is within the interval (-0.6, 2) • lO-^V"^- 

11 MATRE: Most Accurate Time measure- 
ments by REcycled pulsars (Space Watch) 

Exact time measurements are needed in many branches of human activity, 
as well as by simple pedestrians. The nature had prepared us the most 
accurate available watch in the form of Recycled Pulsars. The precision of 
the period of these pulsars is very high, and they are changing their periods 
with time incredibly slowly: dP/dt = IQ-^^-^ for PSR 1953+29, P=6.133 
ms; dP/dt = IQ-^^-^ for PSR 1937+21, P=1.558 ms, what is, for a long 
time period, much better than the highest available precision based on the 
hydrogen frequency standard. 

I suggest to build the international time watch system based on radio 
observations of several most stable recycled pulsars. To build this system 
the following components should be constructed (see fig. 9). 

1. Several radio telescopes over the Earth, which provide stable continu- 
ous observations of the chosen pulsar(s) during 24 hours, and retransmission 
of the signal to the Earth satellites. 

2. System of the Earth satellites (20-25) which accept the exact time 
signals from the radio telescopes, and retransmit them to the Earth, covering 
all the Earth surface, similar to GPS system. 

3 Device which accept the exact tine signal from the satellite, and give 
the exact local/universal time information (should be widely available, could 
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Figure 8: Schematic sketch of the Space Watch system. The signal from 
the recycled pulsar (upper left) is accepted by the Earth radio telescopes, 
which transfer it to the satellite system. The pedestrians accept the exact 
time signal from the satellite by receivers similar to those used in GPS (or 
GLONASS). 
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be combined with the cell phone, or GPS device).). 

4. Software/hardware for transformation of the pulsar signal into the 
signal of the exact time measurement on the radio telescopes, satellites, and 
in every receiving device. 

Such system could be used in all branches of the human activity, including 
ordinary citizen. It should be really international, because its work depends 
on the synchronous and non-interruptive work of all telescopes in different 
countries. The cost of the construction of this system should be quite moder- 
ate (comparable with the construction of GPS), but the exact time is needed 
in much more wide cases than exact coordinates, simply because overwhelm- 
ing majority of people have permanent place of life, or are moving along the 
well known routes. But watches are used by all of them. Therefore, this sys- 
tem could be quite profitable with time, because after construction it needs 
not much for supporting its work. 

12 Conclusions 

1. Timing of the pulsars J0737-3039A /B is the most powerful instrument 
for the verification of General Relativity due to unprecedented precision of 
the observations. 2. Recycled pulsars are the most precise available time 
standards. 3. Checking the physics beyond the standard: G-variability is 
possible by timing of recycled pulsars in close binaries. 4. The service of 
the exact time based on continuous observations of recycled pulsars may be 
constructed and used by pedestrians. 
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